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ABSTRACT 

We report the discovery and follow-up observations of a system of three objects identified by the ALFALFA extragalac- 
tic HI survey, cataloged as (almost) dark extragalactic sources, i.e., extragalactic HI detections with no discernible 
counterpart in publicly available, wide-held, imaging surveys. We have obtained deep optical imaging with WIYN 
pODI and HI synthesis maps with WSRT of the HI1232+20 system. The source with the highest HI flux has a newly 
discovered ultra-low surface brightness (LSB) optical counterpart associated with it, while the other two sources have 
no detected optical counterparts in our images. Our optical observations show that the detected LSB optical counter¬ 
part has a peak surface brightness of ~26.4 mag/arcsec 2 in g ', which is exceptionally faint. This source (AGC 229385) 
has the largest accurately measured HI mass-to-light ratio of an isolated object: Mhi/L g '= 46 Mq/Lq, and has an HI 
mass of 7.2x 10 8 Mq. The other two HI sources (with HI masses 2.Ox 10 8 and 1.2x 10 s M 0 ) without optical counterparts 
have upper limit surface brightnesses of 27.9 and 27.8 mag/arcsec 2 in g ', and lower limits on their gas mass-to-light 
ratio of Mni/Lg/ >57 and >31 Mq/Lq. This system lies relatively close in projection to the Virgo Cluster, but 
velocity flow models indicate that it is located at 25 Mpc, substantially beyond Virgo. The system appears to be quite 
isolated, with no known object closer than ~500 kpc. These HI sources may represent both sides of the threshold 
between “dark” star-less galaxies and galaxies with stellar populations. We discuss a variety of possible formation 
scenarios for the HI1232+20 system. 


1. INTRODUCTION 

Low surface brightness (LSB) galaxies are difficult to 
detect optically, and thus may be underrepresented in 
most optically-selected samples used in studies of galaxy 


history (e.g., Williams et al. 1996 

Madau et al. 1998 

McGaugh et al. 

2000, Brinchmann 

et al. 2004, Hopkins 

Si Beacom 2006 

). Some of the LSB galaxies might be 

those in which si 

,ar formation has been a slow and grad- 

ual process jMcGaugh & de Blok 

1997 

Schombert & 

McGaugh 2(J14p and some may provide a source oi fresh 


tions (Sancisi et al. 1990). By missing the LSB galaxies 
in most surveys, we may be missing an entire population 
of galaxies and/or a phase of galaxy evolution. 

At th e ext reme end of the LSB galaxy spectrum, 
Disney (1976) predicted the existence of entirely “dark 
galaxies , with no observable optical stellar counterparts 
because their surface brightness is too low. A category of 
“crouching giants”, exemplified by the highly luminous 
and massive LSB spiral Malin I (e.g., Lelli et al. (2010)), 
has been identified, but they are quite rare. Overall, no 
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large population of unseen LSB objects has been detected 
at any wavelength. 

LSB galaxies typically possess substantial reservoirs of 
atomic hydrogen, so blind 21-cm surveys represent the 
best opportunity to find large populations of the most 
extreme LSBs. Two major blind HI surveys , HIPASS 
(HI Parkes All Sky Survey, Doyle e t al. |2005b and AL - 
FALFA (Arecibo Legacy Fast ALFA, |Haynes et af.|2l)ll ), 
have reached the conclusion that there is not a significant 
population of gas-bearing but optically dark systems. At 
the same time, there are a number of intriguing, unex¬ 
plained objects detected clearly in HI, showing signs of 
ordered motion and coincident with no discernible stellar 
counterpart The best example of such a “dark galaxy” 
remains the southwestern component of the HI1225 +01 
system ( Chengalur et al.|1995 Matsuoka et al.|2012 |, al¬ 
though it is important to note its presence m a common 
envelope with a visible star forming dwarf companion. 

While some simulations can produce dark galaxies in 
the fo rm of s table gas disks that never produce stars 
(Verde et al. 2002), others find that star-less galaxies 
cannot exist fo r very long before beco ming unstable to 
star formation ( Taylor &; Webster|2005 ). The presence of 
HI in some LSB galaxies provides a ke y dynamical tracer 
of the mass in these extreme systems (Geha et al.||2006" 


Huang et al. 2012aI. Detailed kinematic studies are be- 
ing undertaken to study the effe cts o f outfl ows an d feed - 
back in lower mass galaxies, ( van Eymeren et al.||2009 ), 
in order to understand star formation modes in these 
shallow potential wells and low density galaxies. Groups 
have worked to develop models that can simultaneously 
explain galax y scali ng rela tions in the full cosmological 
context (e.g., jDutton et al~||2007 1. 

Recently, the ALFALFA survey has made significant 
improvements in the sensitivity and depth of available 
wide-held blind HI surveys. ALFALFA has measured 
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Table 1 

Observed HI Parameters of the HI1232+20 System 


AGC 

Position 

Fhi. alfalfa 

Fhi 

V 50 

W 50 

Rhi 

JUxio 19 


J2000 

[Jy km s -1 ] 

[Jy km s -1 ] 

[km s -1 ] 

[km s -1 ] 

['] 

n 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

229385 

12:32:10.3 +20:25:24 

4.87 + 0.04 

4.84 + 0.04 

1348 ± 1 

34+1 

1.60 x 0.72 

1.87 x 1.10 

229384 

12:31:36.4 +20:20:06 

1.36 + 0.03 

1.25 + 0.04 

1309 ± 1 

27+1 

0.73 x 0.55 

1.00 x 0.70 

229383 

12:30:55.3 +20:34:04 

0.81 ± 0.06 

0.42 + 0.05 

1282 ± 4 

59 + 8 


0.80 x 0.28 


Note. — Observed properties of objects in HI1232+20. (1) Catalog ID in the Arecibo General Catalog (an 
internal database maintained by M.P.H. and R.G.) (2) HI centroid position from WSRT (3) Total integrated HI 
line flux density measured from ALFALFA data (4) Total integrated HI line flux density measured from WSRT 
data (5) Heliocentric velocity, measured at the 50% flux level (6) HI velocity width, measured at the 50% flux level 
(7) HI radius at a HI surface density of 1 Mq pc -2 (corresponding to an HI column density of 1.25xl0 20 cm -2 ), 
in arcminutes, measured from the moment 0 maps assuming a beam of 13"x39". Uncertainties on all radius 
measurements are ±0.06'. At a distance of 25 Mpc, 1 ' subtends a distance of 7 kpc. Note that AGC 229384 
is separated into two peaks at the 1 Mq pc -2 level, and that AGC 229383 never reaches a surface density of 
1 Mq pc -2 . Also note that the measurement for AGC 229383 only reflects the radius of the NW clump; the SE 
clump reaches a column density of 5 X 10 1 9 over an area of 14" X 7". (8) HI semimajor and semiminor axis at a 
column density of 5X 10 12 cm -2 , in arcminutes, assuming a beam of 13"x39". 


25,000+ HI sources over 7000 s quare degrees in a cos- 


molo gically si gnificant volume (Giovanelli et al. 2005 
Haynes et al.|20lT Jones et al. in prep.). ALFALFA has 


characterize d thepopulation of normal galaxies (Hu ang 
et al.||2012b|), low mass galaxies (|Huang et al.|20I2a|), as 
well as probing t he HI mass function to lower HI masses 
than ever before (Martin et al. 2010). 

As discussed in Cannon et al. 12015), the ALFALFA 
(Almost) Dark Galaxy Project has been studying the 
very small fraction (^0.4%) of HI sources which lack ob¬ 
vious optical counterparts and are isolated from other 
sources. Followup observations are ongoing and include 
deep optical imaging and HI synthesis maps. Many ob¬ 
jects turn out to be tidal in origin, but some have very 
low surface brightness stellar populations at or below the 
detection limits of current wide field imaging surveys. 

In this work we study the newly discovered HI1232+20 
system of three (almost) dark extragalactic HI sources 
which were not detected in optical surveys, and are at 
least an order of magnitude less lumino us than p reviously 
studi ed LSB galaxy populations (e.g., |Schombert et^al. 
2011). This paper is organized as follows, in Section]5|we 
describe the discovery and observations of this system, 
and in Section [3] we show the results of those observa¬ 
tions. Throughout this work we use a flow model distance 
(Masters 2005) of D =25 Mpc to the HI1232+20 system, 


and we discuss the eff ects of distance uncertainty on our 
conclusions in Section 13.4.21 In Section [4] we discuss the 
implications of these objects and what they might mean 
in the context of (almost) dark galaxies, and in terms 
of extending scaling relationships from normal galaxies. 
Section [5] contains a brief summary of our main results. 
Throughout this work we assume a ACDM cosmology, 
with fl m =0.3, Da=0.7, and H o =70 km s -1 Mpc -1 . 

2. OBSERVATIONS 

2.1. ALFALFA discovery of HI1232+20 system 

The ALFALFA survey employs a two-pass, fixed az¬ 
imuth drift scan strategy, the details of which are de¬ 
scribed in previous paper s (Giovane lli et al.||2 005 Sain- 
tonge|[2007 Martin et al.||20U5 ' Haynes 5 ^Fm!pOTl ). All 


data are flagged tor radio frequency interference (RFI) 
interactively, and each grid is examined by hand to con¬ 
firm and improve on sources detected via the automated 


methods of Saintonge |2007); final source parameters are 
measured and cataloged interactively. 

Among the ALFALFA (almost) dark extragalactic 
sources, the HI1232+20 system (comprised of sources 
AGC 229383, AGC 229384, and AGC 229385) was found 
to be of particular interest. These three objects are near 
each other on the sky and also have similar recession ve¬ 
locities, so are likely associated with each other. From 
the ALFALFA observations it was clear that these three 
sources have significant amounts of gas present, even 
though they do not have readily identifiable optical (stel¬ 
lar) counterparts in existing optical databases (SDSS, 
DSS). While they appear on the sky near AGC 222741 
(CGCG 129-006), there is a significant separation in ve¬ 
locity between the sources. AGC 222741 has an HI re¬ 
cession velocity of 1884 km/s while the three sources in 
this sample have recession velocities of ~1300 km/s. 

An overlapping archival ultraviolet (UV ) image fro m 


GALEX GR7 (Galaxy Evo lution Explorer, Martin et al. 
™ |Morrissey et al. 2007 , Data Release 7, |Bianchi et al. 


2005 


2014) shows a faint diffuse UV source at the coordinates 


of ACC 229385 (see Section 2.4 for more details). There 
is also a hint of a faint object at the same position in 
the DSS2-B image (Digitized Sky Survej®| but no source 
visible at that position in images from SDSS DR9 (Sloan 
Digital _Sky Sjurvey, Data Release 9, Ahn et al. 2012). 
Tables [T| and [2] contain complete information about the 
HI1232+20 system. No optical sources were evident at 
the locations of the other two HI detections. 

Given the curious nature of this system, we have car¬ 
ried out further observations to study it in more detail. 
We have obtained deep optical images to look for possi¬ 
ble faint stellar populations in the sources, and sensitive 
HI synthesis observations to resolve the gas distribution 
and kinematics in more detail. 


8 The Digitized Sky Surveys were produced at the Space Tele¬ 
scope Science Institute under U.S. Government grant NAG W- 
2166. The images of these surveys are based on photographic data 
obtained using the Oschin Schmidt Telescope on Palomar Moun¬ 
tain and the UK Schmidt Telescope. The plates were processed 
into the present compressed digital form with the permission of 
these institutions. 
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Figure 1. Detection image made from all pODI observations with ALFALFA HI sources in the HI1232+20 system labeled and WSRT 
contours overlaid in red. The size of the ALFALFA and WSRT beams are indicated by the shapes at the top left corner. The WSRT 
pointings cover nearly the entire area region shown in this image. The optical image is binned 12 X 12 pixels so that the resulting pixels 
are 1.3” on a side, and is aligned north up east left. The WSRT contours show column densities spaced logarithmically between 1 and 
64xl0 19 cm -2 . The tight packing of the lowest contours is likely a result of sigma-clipping, and not a sharp edge in the HI distribution. 
See Section |2.3| for more details and about sensitivity differences across the image. The spiral galaxy AGC 222741 has a recession 
velocity of cz= 1884 km/s, s o is sig nificantly more distant than the objects in the HI1232+20 system, which have recession velocities around 
cz~1300 km/s (see Section |3. 4. 2| for details). Th e W SRT contours of AGC 229383 show that it is separated into two components at this 
sensitivity level, as is discussed further in Section |3-3| Also visible at the top left of the image is diffuse filamentary emission from Galactic 
cirrus, which coincides with features seen in far-infrared and ultraviolet images of the same area. 


2.2. Deep optical imaging with WIYN pODI 

The HI1232+20 system was observed with the WIY1N0 
3.5-m telescope at Kitt Peak National ObservatorjFj us¬ 
ing the partially populated One Degree Imager (pODI). 
Currently, pODI is made up of 13 Orthogonal Trans- 

9 The WIYN Observatory is a joint facility of the University of 
Wisconsin-Madison, Indiana University, the University of Missouri, 
and the National Optical Astronomy Observatory. 

10 Kitt Peak National Observatory, National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy (AURA) under cooperative agreement 
with the National Science Foundation. 


fer Arrays (OTAs), each of which is made of sixty-four 
480x496 pixel cells. The OTAs are arranged on the fo¬ 
cal plane such that the central 3x3 OTAs cover an area 
of 24' x 24' with pixels that are 0.11" on a side. The 
numerous gaps between cells and OTAs require a series 
of offset dithered exposures to produce a well-sampled 
image. Four of the standard Sloan Digita l Sky Survey 
(SDSS,|Gunn et al.l (|1998[) , |Doi et al.l 12010b) g',r',i', and 
ab 


z' filters are available, and stars from the SDSS catalog 
photometry are used for standard photometric calibra¬ 
tions. 

We imaged an area which includes both AGC 229384 
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and AGC 229385 on the night of 6 February 2013 with 
nine dithered 300 second exposures in each of the g'r'i' 
filters. We observed AGC 229383 on 2 May 2014 with 
nine dithered 300 second exposures in both the g' and 
r' filters. By combining data from these two nights of 
observations, we have contiguous deep multi-wavelength 
imaging coverage over an area ~40' x 40'. We also im¬ 
aged this field with an 80A narrow-band Ha filter during 
photometric conditions on the night of 6 February 2013. 
Our dithered sequence of nine 300s images, while not cal¬ 
ibrated, do not show any Ha detections at the locations 
of the three HI sources in the HI1232+20 system, but 
do show a background spiral galaxy (AGC 222741) quite 
clearly. 

We reduced our observations using the QuickReduce 
(QR, Kotulla| 2014) data reduction pipeline, and supple- 
mented this processing with an additional illumination 
correction. QR was run interactively in the One Deg ree 
Imager Pipeli ne, Portal, and Archive (ODI-PPA W Tsci- 
ence gateway ( Gopu et al.|2014 Young et al.|2 013'b The 
PPA interface allows the user to select which observa- 
tions will be reduced, and runs all of the reductions on 
computing resources at the Pervasive Technology Insti¬ 
tute (PTI) at Indiana University. 

The QR pipeline includes: masking of saturated pix¬ 
els, crosstalk, and persistence; overscan subtraction; bias 
level subtraction; dark current subtraction; non-linearity 
corrections to each cell; flat Held correction from dome 
flat fields; cosmic-ray removal; fringe removal (in *'); 
pupil ghost correction. However, the final pipeline- 
processed data still have uncorrected instrumental ar¬ 
tifacts in them, especially at very faint intensity levels. 
In order to produce images that are suitable for low sur¬ 
face brightness analysis, we need to correct for the small 
gradients, sky level offsets, and other artifacts in partic¬ 
ular cells. Once these effects are corrected, the dithered 
images can be combined into a final deep image. 

In order to remove these image artifacts, we apply an il¬ 
lumination correction using dark sky flats generated from 
the observations themselves. For a particular filter, we 
mask all objects in the images, then use a median algo¬ 
rithm to combine all of the exposures into a dark sky flat 
field, which is then smoothed with a 3 x 3-pixel smoothing 
element. Each exposure is then divided by this illumina¬ 
tion correction image. 

Before combining all exposures in a dither pattern, we 
re-project them to a common pixel scale and also scale 
the images to a common flux level using measurements 
of stars in t he fie ld and SDSS DR9 catalog magnitudes 
(Ahn et al. 2012). This compensates for varying sky 
transparency during the dither sequence, and typically 
yields final photometric zeropoints with standard devia¬ 
tions of 0.02 — 0.03 magnitudes. The g' filter calibrations 
required a g'-i' color term of amplitude 0.079±0.013, but 
r' and i' calibrations required no color term. The point 
sources in our final combined images have an average 
FWHM of 0.7" - 0.9". 

We also create a deep “detection-only” image by com- 


11 The ODI Pipeline, Portal, and Archive (ODI-PPA) is a joint 
development project of the WIYN Consortium, Inc., in partnership 
with Indiana University’s Pervasive Technology Institute (PTI) and 
with the National Optical Astronomy Observatory Science Data 
Management (NOAO SDM) Program. 


bining all images from both pointings in all filters to 
reach the faintest light possible. This detection image 
is then binned to 1.2" resolution to bring out very faint 
emission, and is shown in Figure [I] with relevant HI 
sources labeled and HI synthesis contours overlaid. The 
contrast levels in this image have been stretched to show 
the exquisite sensitivity to faint light. In this view, the 
optical counterpart to AGC 229385 is strikingly visible, 
as will later be discussed. Also visible in the upper left 
corner of the image is diffuse filamentary emission from 
Galactic cirrus. This foreground emission comes from re- 


(Sandage 

1976,[Witt et al. 

2008). 

vevs (IRi 5 

lS, Schiegei et al 

1998 


show features that are coincident with the faint optical 
emission we see in our image. It even shows up weakly in 
the overlapping deep archival GALEX UV image. Galac¬ 
tic cirrus is very faint and diffuse at optical wavelengths 
and typically only visible in deep, w ide-field images tha t 
are very accurately flat-fielded (e.g., iRudick et al.|2010). 


2.3. HI synthesis imaging with WSRT 

We observed the HI1232+20 system with four 12h 
pointings at the Westerbork Synthesis Radio Tele¬ 
scope (WSRT), three of which were centered on 
(12:31:52.6 +20:22:59) to encompass the centroids of 
AGC 229385 and AGC229384, and one of which was 
centered on (12:31:08.5 +20:31:41.9), to encompass 
AGC 229383. The primary beams of the two pointing 
centers are 35' wide and cover nearly all of the area dis¬ 
played in Figure [lj We observed the HI line in one band 
with 10 MHz bandwidth, two polarization products, and 
1024 channels, ensuring a broad range of line free chan¬ 
nels for continuum subtraction and a velocity resolution 
of 4.12 km s -1 after Hanning smoothing. 

The data were reduced using the same aut omated data 
reduction pipeline as applied in Wang et al. (2013), orig¬ 


inally used by Serra e t al. 
tion software Miriad (|Sault 


(2012t 

t et al. 


using the data reduc- 


19951 wrapped into a 


Python script. The data were automatically flagged for 
radio interference using a clipping method after filtering 
the data in both the frequency- and time-domain. Af¬ 
ter the primary bandpass calibration, the data were it¬ 
eratively deconvolved with the CLEAN algorithm, using 
clean masks determined on the cube with decreasing clip 
levels, to then apply a self-calibration. The calibration 
solution was applied to the visibilities and the continuum 
was subtracted in the visibility domain to then invert the 
data after Hanning smoothing, using a set of combina¬ 
tions of Robust weighting and tapering with a Gaussian 
kernel, as well as a binning in the frequency domain. Fi¬ 
nally the data cubes were iteratively cleaned using clean 
masks determined by filtering the data cubes with Gaus¬ 
sian kernels and applying a clip level. The clean cutoff 
level was set to the rms noise in the data cubes. Because 
we hence cleaned the data comparably deeply, no correc¬ 
tion of the intensity levels of the residuals was made. 

Our pipeline produces cubes at each centroid with 
three different robustness weightings, r = 0.0, r = 0.4, 
and r = 6.0, binned to a velocity resolution of 6.2 
km s -1 after Hanning smoothing (12.4 km s -1 for the 
r = 6.0 cube). The noise level in the cubes for the 
three respective robustness weightings are 0.40, 0.36, and 
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Figure 2. Comparison of WIYN pODI image and GALEX FUV image. The pODI detection image is binned 12 X 12 so that pixels are 
1.3" on a side. The FUV image has been binned 3 X 3 so that pixels are 13” on a side. Two of the three HI sources in the HI1232+20 
system are labeled on both panels with 3 7 circles, the size of the ALFALFA beam. AGC 229385 clearly has an optical and ultraviolet 
counterpart while AGC 229384 lacks a counterpart at either wavelength. AGC 229383 has not been observed by GALEX. 


0.24 mJy/beam/channel, with beam sizes of 39" xl3", 
45"xl5", and 54"x20". 

For each cube we then created HI total flux maps 
by summing masked cubes along the velocity axis. We 
created the masks by smoothing the images to twice 
the beam size, and then keeping any pixel 3 a above 
the noise level. From these we calculate HI column 
density maps assuming optically thin HI gas such that 
N#/=1.823 x 10 18 / Tbdv cm -2 . Since the final contour 
map results from the combination of multiple WSRT ob¬ 
servations (3 at the SE pointing, and 1 at the NW point¬ 
ing), the signal-to-noise ratio varies across the image and 
is less sensitive near AGC 229383. The lowest HI contour 
shown on the images, lxlO 19 cm -2 , corresponds to a less 
significant detection in the region around AGC 229383 
than it does in the region around AGC 229385 and 
AGC 229384. As a result, there were locations outside of 
the main signal from AGC 229383 where the HI column 
density exceeded lxlO 19 cm~ 2 , but since the significance 
of the detection was <3er, those contours are not shown. 

We additionally create a one dimensional integrated HI 
line profile for each object, as displayed in Section [3] We 
fitted the line with both the two horned function applied 
in the ALFALFA data processing, and using a standard 
Gaussian fit and note that the fluxes from the fits match 
well within random errors. We recover 99% of the AL¬ 
FALFA flux in the WSRT spectrum of AGC 229385, and 
92% in AGC 229384, but only 52% in AGC 229383. The 
spectra for AGC 229385 and AGC 229384 are both well 
fitted by a Gaussian profile, and though both may show 
slight deviation from Gaussian, using two Gaussians does 
not return a better result. The spectrum for AGC 229383 
is not well fitted by either a Gaussian or a two-horned 
fit. 

Finally, we produce HI velocity maps using two differ¬ 
ent methods. We created standard moment 1 maps from 
cubes masked at 3cr, and additionally fitted Gaussian 
functions to each individual profile in the datacube us¬ 
ing the GIPSY task XGAUFIT. The resulting maps from 


the two methods are virtually identical, and we show the 
maps in Section [3] 

To further analyze the velocity field we create Position- 
Velocity (PV) diagrams for each source. We produced 
the PV diagrams by taking an 18" (1 minor axis beam 
width) wide slice along the HI major axis centered on the 
HI surface density centroid. We measured the position 
angle and centroid from the surface density profile since 
the variations in the velocity field leave the major veloc¬ 
ity axis and center uncertain. We note that none of the 
PV diagrams change significantly for small variations in 
position angle or slice width. 


2.4. Archival GALEX observations 


AGC 229385 (the strongest HI detection of the 
HI1232+20 system) had a very faint UV co unterpart vis¬ 


ible in an archival dataset from GALEX (Martin et al., 
2005 Morrissey et al.|2007). GALEX obtained images in' 


the far ultra-violet (FUV) from 1344 to 1786A with 4.3" 
FWHM resolution, and in the near ultra-violet (NUV) 
from 1771 to 283lA with 5.3" resolution. These UV im¬ 
ages are especially sensitive to young stellar populations, 
and should help to identify sites of recent star formation. 

AGC 229385 was imaged by GALEX in May 2007 
in the NUV and FUV bands, with exposure times of 
1145s in both bands. A set of matched images is shown 
in Figure [2j with data from pODI WIYN and GALEX 
FUV. The locations of AGC 229385 and AGC 229384 


are shown in these images, while AGC 229383 lies out¬ 
side of any archival GALEX image. The brightest HI 
source, AGC 229385, is faintly visible in the UV image 
as a diffuse source. However, the GALEX pipeline (GR7, 
Bianchi et al. 2014) does not identify this diffuse object as 
a source, instead shredding it into multiple point sources. 
We measure the brightness of AGC 229385 in the FUV 
and NUV images ourselves in an aperture matched to 
our optical images, as discussed in Section [3] There is 
no source visible in the FUV image at the position of 
AGC 229384. 
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Table 2 

Derived properties and limits from observations 


Quantity [units] 

AGC 229383 

AGC 229384 

AGC 229385 

m g , [mag] 

> 20.7 

> 20.7 

19.20 (0.03) 

m r r [mag] 

> 20.2 

> 20.5 

19.27 (0.03) 

[mag] 


> 19.7 

19.36 (0.04) 

Mg',peak [mag/arcsec 2 ] 

> 27.8 

> 27.9 

26.4 (0.1) 

Mr', peak [mag/arcsec 2 ] 

> 27.3 

> 27.7 

26.5 (0.1) 

Mi', peak [mag/arcsec 2 ] 


> 26.8 

26.1 (0.1) 

Optical major axes (at fJ, g / = 27) [kpc] 



7x3 

HI major axes (at 5xlO l9 cm 2 ) [kpc] 

12 x 4 

14 x 10 

28 x 16 

M g , [mag] 

> -11.3 

> -11.3 

-12.89 

M r , [mag] 

> -11.8 

> -11.5 

-12.79 

M v [mag] 


> -12.3 

-12.69 

g' - r' [mag] 



-0.09 

B — V [mag] 



0.13 

M b [mag] 

> -11.2 

> -11.1 

-12.72 

Lfuv [£©] 



2.41 x TO 7 

Lnuv [£©] 



2.66 x 10 7 

SFR N uv [M©/yea,r] 



4.1 x 10“ 3 

SFR FU v [M©/year] 



6.9 x 10“ 3 

log M H i [log Mq\ 

8.08 

8.30 

8.86 

M* [M 0 ] 

<3.7 x 10 5 

<3.4 x 10 s 

1.5 x 10 6 

Mh 1 /Mi, 

> 320 

> 580 

290 

M HI /L g , [Mq/Lq] 

> 31 

> 57 

45.8 

Mhi/Fb [Mq/Lq\ 

> 26 

> 48 

38.2 

Note. — Apparent magnitudes (m g < 

mg) 

are not corrected 

for Galactic 


extinction. Absolute magnitudes ( M g /, M r /, M^/, Mb) luminosities (LpxjVi L nuv)i 
and colors ( q' — r\ B — V) are corrected for Galactic extinction from Schlafly &| 
|Finkbeiner|(|2011|). All absolute quan tities assume a dist ance of 25 Mpc. Mb and B — 
V are determined from conversions i n|Jester et al. ]( |2005| . Upper limits are determined 
where sources are not detected in pOJJl observations and are at 3a confidence levels. 
Uncertainties on measured quantities are indicated in parentheses. 

3. RESULTS 

In the following sub-sections we describe the results 
of our followup observations for the HI1232+20 system. 

The derived results are summarized in Tabled We also 
consider the environment around this system, and uncer¬ 
tainties in the adopted distance. 

3.1. AGC 229385 

Figure [3] shows our observations of AGC 229385. Fig¬ 
ure [3ja is a deep 3-color image from WIYN pODI with 
contours from our WSRT HI synthesis map at the highest 
resolution (39"xl3"). Figure [sf shows a zoomed-in re¬ 
gion around the optical counterpart with the regions we 
use to measure its surface brightness indicated by black 
squares. The optical emission from AGC 229385 appears 
very blue, and coincides spatially with the peak of the HI 
distribution. The optical component appears ~5x less 
extended than the radio emission, but both are similarly 
elongated in the northeast-southwest direction. 

AGC 229385 has an unusual optical morphology that is 
not simple to describe. The optical counterpart is elon¬ 
gated in the N-S direction, but has a nearly constant 
surface brightness across its entire extent. We fitted an 
ellipse to the 5 a contour on the g' image and found a 
semi-major axis of 32" and semi-minor axis of 10". This 
5er ellipse has a position angle of 15° (measured clockwise 
from N) and an ellipticity (e=l — b/a) of e=0.68. 

After attempts to fit elliptical annuli to the optical 
images of AGC 229385 resulted in inconsistent and di¬ 
vergent surface brightness profiles, we decided to mea¬ 
sure the surface brightness in small regions instead. We 
measured the optical surface brightness of AGC 229385 


in 5" x 5" regions following the curving shape of this 
source from south to north, after masking all obvious 
foreground and background sources. These regions are 
shown as black boxes in Figure [3Jd. We also placed simi¬ 
lar regions outside the periphery of the source to deter¬ 
mine the local sky value. The surface brightness traces 
along AGC 229385 are plotted in the Figure [3}© We cal¬ 
culated the surface brightness level that corresponds to 
3 times the standard deviation in the sky, and label it 
as the 3 a detection threshold. The g' surface brightness 
profile has the highest signal-to-noise ratio, and is well 
above the 3er level. The r' profile is weaker but still well- 
measured. In i', this source is only weakly detected, but 
still is above the 3cr surface brightness level. In all three 
filters a similar profile shape is seen as the outlying re¬ 
gions show very little signal and the inner regions show 
a relatively flat brightness distribution across the source. 
To estimate a peak surface brightness value in each fil¬ 
ter we use the measurements from the three boxes just 
south of the bright foreground star near the center of 
AGC 229385, where there are relatively few contaminat¬ 
ing sources which had to be masked and the profiles are 
relatively smooth. The peak values are calculated by av¬ 
eraging the measurements in these three boxes, and are 
26.4, 26.5, and 26.1, in g',r', and i', respectively. While 
the formal uncertainties on these surface brightness mea¬ 
surements are low (~3 — 5%, owing to our accurate pho¬ 
tometric calibrations and the good S/N of the optical 
counterpart), the variations between adjacent boxes can 
be as high as 0.1 mag. Accordingly, we assign an uncer¬ 
tainty of 0.1 mag to these peak values. We also measure 
photometry of this source in a 32" radius aperture af- 
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Figure 3. Our observations of AGC 229385. a) 3-color optical image from WIYN pODI binned 4 x 4 so that pixels are 0.44" on a 
side. North is up, East is left, and the angular scale is indicated on the image. White contours of HI emission are from WSRT and are 
at are at 1, 2, 4, 8, 16, 32, and 64 x 10 19 atoms cm -2 , and the white ellipse at the bottom left shows the size and shape of the WSRT 
beam, b) Zoomed-in view of the WIYN pODI observations, with black boxes indicating regions used for surface brightness measurements. 

c) Surface photometry for AGC 229385 from pODI images in g ', r', and i' filters, measured in b n x 5" boxes following the brightest part 
of the optical emission, after masking obvious stars and background galaxies. Dashed lines show the 3cr surface brightness detection limits 
in each filter. Blue points and lines indicate g' filter measurements, green indicate r' measurements, and red indicate i' measurements. 

d) XGAUFIT velocity field from WSRT with the same HI contours as above. The ellipse in the bottom left corner indicates the size and 
shape of the WSRT beam. The dashed rectangular region shows the slice along the major axis which is used to generate the P-V diagram. 

e) Position-Velocity diagram from WSRT. f) ALFALFA and WSRT HI spectra of AGC 229385. As discussed in the text, HI fluxes from 
ALFALFA and WSRT are in good agreement. Blue dashed lines represent the ALFALFA spectra, while the solid black lines show the 
WSRT r = 0.4 spectra. 


ter masking obvious foreground and background sources. 
The results of the surface and aperture photometry of all 
three sources are summarized in Tabled] 

The HI synthesis observations from WSRT (shown in 
the bottom row of Figure [3]) are also difficult to inter¬ 
pret. The moment 0 map, Flgure[3]l, shows an HI source 
significantly more extended than its optical counterpart. 
AGC 229385 has HI major and minor axes of 3.2'x 1.4' 
measured at an HI surface density of 1 M 0 /pc 2 (corre¬ 
sponding to a column density of 12.5 x 10 19 cm 2 ), or 
24 kpc x 10 kpc assuming a distance of 25 Mpc. WSRT 
measures significant emission at lower column densities, 
out to 3.7' x 2.2' or 28 kpc xl6 kpc at 5 x 10 19 cm -2 and 
a furthest extent of ~5' at 10 19 cm -2 . 

The exact surface density profile of AGC 229385 is sub¬ 
ject to its 3D geometry and inclination, which are diffi¬ 
cult to determine conclusively given our beam size and 
the unusual nature of AGC 229385. AGC 229385 ap¬ 
pears to have an HI position angle of 21° which gives 
~6 x 7 resolution elements at the furthest extent along 
the major and minor axes. We formally measure an in¬ 


clination of 63 ± 4 degrees assuming a thin HI disk and 
uncertainties of half the beam size along the major and 
minor axes. 


As an instructive exercise, we assume a disk geom¬ 
etry and compute de-projected surface density profiles 
for AGC 229385 using Ro bertson-Lucy deconvolution 


(jLucy 

1974 

Warmels 

1988 

the GIPSY task RADIAL) 

This method 

, developed tor use in low resolution imag- 


ing, works by collapsing the measured intensity along 
the minor axis, and produces a one-dimensional pro¬ 
file, which is then iteratively matched by a model one- 
dimensional profile produced by summing axisymmetric, 
uniform density co-planar rings along lines of sight. This 
method does not require knowledge of the inclination 
of the object, but still assumes a disk geometry. The 
summed one-dimensional surface density profile shows 
some asymmetry and two peaks with a slight depression 
in the center. These features are reflected by asymmetry 
in the resulting RADIAL model, and a strong suggestion 
of a hole in the center of the HI distribution, a feature 
that would be smeared out in 2D ellipse fitting analysis. 
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Figure 4. Our observations of AGC 229384, same panels as Figure 13] No optical counterpart is detected down to an upper limit of 
fi g i ~27.8 mag/arcsec 2 , so no surface brightness profiles are shown. The black “x”s in panels a and b mark the two peaks of the HI 
distribution. The HI velocity field of AGC 229384 in panel e shows no evidence of ordered rotation. As with AGC 229385, the WSRT and 
HI spectra in panel f are in good agreement with each other. 


If confirmed, this hole could be indicative of the forma¬ 
tion of cold atomic and molecular hydrogen in the center 
of the object, or of a non-disky, more complicated HI dis¬ 
tribution, possibly caused by two recently merged com¬ 
ponents. However, higher resolution data are necessary 
to confirm the existence of a gap in the HI distribution. 

The HI velocity held of AGC 229385 shown in Fig¬ 
ure w is equally difficult to interpret. The narrow inte¬ 
grated HI line width and single peaked spectrum (shown 
in Figure J3t) is suggestive of slow rotation. If we assume 
a thermalvelocity dispersion of 11.0 km/s, the 34 km/s 
integrated line width of AGC 229385 gives an observed 
rotation velocity of 16 km/s when subtracting thermal 
velocity and dividing by 2. This translates to a rough 
inclination-corrected rotation velocity of 18 km/s. In¬ 
deed, the moment 1 map shown in Figure [3j:l shows evi¬ 
dence of ordered rotation roughly along the major axis of 
the source, but the gradient is asymmetric/The irregular 
shape of the P-V diagram for AGC 229385 (Figure^) 
further diagnoses this asymmetry. The southern side of 
the source shows a clear slope of 12 — 15 km/s, but then 


any turnover in the north side of the object, but does not 
give a significant gradient. 

The FUV image of AGC 229385 from GALEX is shown 
in Figure [2] While this image is less striking than the 
optical images, this source is still detected in both the 
NUV and FUV images. The GALEX pipeline shreds 
this extended source into multiple point sources, so we 
measure its brightness in the same 32" radius aperture 
and with the same masking that was used on the op¬ 
tical images. The apparent magnitude and uncertain¬ 
ties in the NUV and FUV bands are 19.631 (0.069) mag 
and 19.155 (0.035) mag, respectively. After correcting 
for Galactic extinction, we use the assumed distance of 
25 Mpc to determine star formation rates (SFRs) from 
the NUV and FUV luminosi ties following th e relations of 
Murphy et al.l (|2011h and |Hao et al.| (|201 1|) , and report 
the results in Table (5J We ao~nob include a correction for 
internal extinction. 

Using the flow model distance of 25 Mpc, we derive 
absolute global parameters for AGC 229385, which are 
listed in Table [2| AGC 229385 has an optical luminos- 


one approaches the north side of the galaxy. Further, 

Mb=— 12.72 via Jester et al. 

2005 

). Its g' - 

- r' color 

the P-V diagram reveals that the velocity dispersion is 

is very blue, and corresponds 

to a 13 — V color o± U.13 

of a similar order to the velocity gradient. It is possible 

(via similar conversion in pester et al. 

2005 

). Using 


density major axis: fitting a P-V diagram at 0° removes we find Mhi/Lb=38.2Mq/Lq. Using the self-consistent 
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Velocity (km/s) 


Figure 5. Our observations of AGC 229383, same as Figure [3] a) WIYN pODI optical image covering both peaks of the HI distribution, 
with contours from WSRT and black “x”s marking the HI peaks, b) A zoomed in view of the NW component, with a black “x” marking 
the peak of the HI distribution. No optical counterpart is detected. An asteroid trail is also visible in this image, c) A zoomed in view 
of the unresolved SE component, with its center marked. Again no optical counterpart is detected, d) The moment 1 map from WSRT 
with the same HI contours as above. The ellipse in the bottom right corner indicates the size and shape of the WSRT beam. The dashed 
rectangular region shows the slice along the major axis which is used to generate the P-V diagram to the right, e) Position-Velocity diagram 
from WSRT. No ordered rotation is visible in this splotchy and irregular velocity field, f) ALFALFA, WSRT, and LBW HI spectra of 
AGC 229383. As discussed in the text, LBW recovers 64% of the ALFALFA HI flux, and WSRT recovers only 52%. This likely means 
that there may be extended gas below the WSRT sensitivity which connects these sources, and we consider them to be peaks of a common 
source in this work. 


simple stellar popula tion models discussed 


in |McGaugh] 
stellar mass- 


& Schombert (2014), this color implies a 
to-light ratio (in the V band) of M*/Ly=O.3M 0 /L 0 , 
or a total stellar mass of ~1.5 x 10 6 M 0 , and a ratio 
f H i=M HI /M+= 475. 


3.2. AGC 229384 

Figure |4ji shows the deep 3-color WIYN pODI image 
of AGC 229384 with WSRT HI contours overlaid. No 
obvious optical counterpart is visible. The faint grid of 
horizontal and vertical stripes in the background is an ar¬ 
tifact from the data reduction process. Figure [4 )d shows 
a zoomed-in view near the twin peaks of the HI con¬ 
tours (marked with black xs). We used small 5" x 5" 
regions placed around this area to determine upper lim¬ 
its on the optical non-detection, using the same method 
as for AGC 229385, after masking all obvious foreground 
and background sources. The 3er upper limits on this 
non-detection are given in Table [2j and are 27.9, 27.7, 
and 26.8 mag/arcsec 2 in g\ r', and i', respectively. To 
estimate an upper limit on an integrated magnitude, 
we must assume an aperture size. Since the HI major 
axis of AGC 229384 is ~50% of the HI major axis of 
AGC 229385, the aperture is scaled by the same amount, 
to 15". Using this aperture of radius 15", we find 3 a 


upper limits on the integrated magnitude in g ', r', and 
i’ filters to be 20.7, 20.5, and 19.7 mag, respectively. 
These upper limits are used to generate upper limits on 
Mhi/Lq 1 , M*, and Mhi/M *, all of which are shown in 
Table [2] In particular, these non-detections correspond 
to a stellar mass upper limit of M*<3.4 x 1O 5 M 0 . When 
necessary, the observed colors of AGC 229385 were used 
to make filter conversions on the optical non-detections 
for AGC 229384. An archival GALEX image that covers 
AGC 229384 also shows no optical counterpart for this 
object (see Figure [2]). 

The HI contours of AGC 229384 in Figure [4jr show 
an irregular distribution with two weak density peaks. 
The HI major axes at a surface density of 1M 0 pc -2 are 
1.4' x 1.2' or 10 kpc x 8 kpc. The HI velocity field shown 
in Figure |4ji is patchy and irregular, and appears to be 
dominateaby random motions. The position-velocity di¬ 
agram in Figure® shows no evidence of ordered rotation. 
At its assumed distance of 25 Mpc, AGC 229384 has a 
total HI mass of Mhi= 2.0 x 1O 8 M 0 . Complete details 
are given in Tables |T] and [2j 

3.3. AGC 229383 

From the ALFALFA HI observations, AGC 229383 was 
extracted as a single, possibly extended weak source. Fol- 
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low up observations with the single pixel L-Band Wide 
(LBW) receiver confirmed the existence and extended 
nature of the source, since the more sensitive LBW ob¬ 
servations only recovered 64% of the original ALFALFA 
flux, as shown in Figure [5j\ WSRT observations resolve 
the source into two low HI column density clumps, sep¬ 
arated by 5.5 arcminutes (~ 40 kpc at D = 25 Mpc). 
The SE clump is only detected in a single beam at low 
signal to noise ratio, but is detected in both WSRT point¬ 
ings which overlap its position. It is possible that these 
two clumps are independent sources. However, the two 
clumps together only recover 52% of the original AL¬ 
FALFA flux, suggesting that there may be gas connecting 
the sources below the sensitivity of the synthesis obser¬ 
vations. AGC 229383 may be two distinct sources, but 
because the missing flux and the clumpy HI distribution 
are ambiguous as to the true nature of the source, we 
choose to discuss it as a single source with two peaks in 
the remainder of this paper. 

Figure [5^, shows the color image made from the g' and 
r' observations of AGC 229383, which is located 19' to 
the NW of AGC 229385. As with AGC 229384, no opti¬ 
cal counterpart is visible. Figure [5b shows the zoomed- 
in region around the NW peak of the HI distribution. 
Again we masked obvious foreground and background 
sources and used regions around this area to measure 
the background statistics and determine upper limits on 
the optical non-detection using the same method as for 
AGC 229385. These 3er upper limits are given in Ta¬ 
ble [ 2 ] and are 27.8 and 27.3 mag/arcsec 2 , in g' and r', 
respectively. Similarly, we find upper limits on integrated 
magnitudes in a 15" aperture of 20.7 and 20.2 mag in 
g' and r', respectively. We use these in the same way 
as AGC 229384 to generate upper limits for the derived 
quantities for AGC 229383, all of which are listed in Ta¬ 
ble^ Figure [5J; shows the zoomed-in region around the 
unresolved SE peak of AGC 229383, and reveals no op¬ 
tical counterpart as well. 

The HI contours of AGC 229383 in Figure [5^t show the 
two density peaks, of which the NW component appears 
resolved and the SE component unresolved. The HI 
never reaches a projected HI surface density of lM^pc -2 
(corresponding to a column density of 12.5 x 10 19 cm -2 
at this distance), and has a maximum extent of only 
48" x 17" at a column density of 5xl0 19 cm -2 . The HI 
velocity field is patchy and irregular, and appears to be 
dominated by random motions, though our data are lim¬ 
ited due to the fact that the two peaks are poorly resolved 
or unresolved. The limited data and patchy nature of the 
velocity field precluded meaningful fitting with XGAU- 
FIT, so instead we show the moment 1 map in Figure[5tl, 
and the messy position-velocity diagram in Figure [5j\^t 
its assumed distance of 25 Mpc, AGC 229383 has a total 
HI mass of 1.2 x 10 8 Mq. 


3.4. Isolation, Environment, and Distance Uncertainty 
3.4.1. Isolation of the HI1232+20 system. 

Since many “dark” galaxy candi dates turn out to 
be tid al features (e.g., VIRGOHI21 |Duc & Bournaud 
(2008)) rather than isolated galaxies, we look tor pos- 
sible objects which may have recently tidally interacted 
with HI1232+20. To test this possibility, we searched all 
cataloged nearby sources in the Arecibo General Cata- 



Figure 6. All panels show galaxies from the UZC, SDSS, and 
AGC. Larger circles represent more luminous galaxies, and smaller 
circles and dots represent less luminous galaxies. Note that we 
are plotting the three objects in the HI1232+20 system as red 
dots, as well as the background spiral galaxy (AGC 222741) at 
C2:=1884 km/s as a red “x”. The latter is nearby the objects in 
the HI1232+20 system in the RA-Dec plane, but is well-separated 
in the other two panels. The bottom left panel shows the galaxy 
distribution in the plane of the sky. The top left panel shows 
the galaxy distribution in the velocity-right ascension plane. The 
bottom right panel shows the galaxy distribution in the declination- 
velocity plane. The solid curved lines show the virial radius of the 
Virgo Cluster in each projection. Dotted circles enclose the galaxies 
nearest to this system. 

log (AGC), NED, and the SDSS spectroscopic survey and 
determined the timescale on which they could have in¬ 
teracted with this system, given their current velocities. 
If we generously assume that a flyby encounter may have 
had a relative velocity of ~500 km/s, we can calculate 
how long ago the nearby objects could have interacted. 
Naturally, the three sources in the HI1232+20 system all 
have short interaction time scales with each other (<300 
Myr, based on these assumptions). The only source with 
an interaction time scale <1 Gyr is AGC 742390, which 
is ~30' W of this system and has cz = 1127 km/s. 
AGC 742390 is an elongated star-forming galaxy with 
M g = —15 mag, assuming it is at the same distance of 
25 Mpc. At this distance, AGC 742390 has a projected 
separation from the HI1232+20 system of ~200 kpc, but 
its smaller recession velocity (Av ~ 200 km/s) implies 
that it is likely more nearby than this system. Given the 
lack of obvious nearby galaxies in optical and HI surveys, 
and the lack of objects which could have recently tidally 
interacted with this system, it seems that HI1232+20 is a 
locally isolated system, and not a tidal feature of a larger 
parent object. Still, we cannot exclude the possibility 
that this system may have been produced as a result of 
tidal interactions or other gas stripping processes. 


3.4.2. Effects of distance uncertainty and the environment 

Througho ut this wor k we have adopted a flow-model 
distance (JMasters|2005 ) of D =25 Mpc to the HI1232+20 
system. However, given its location on the outskirts of 
the Virgo Cluster, there is some uncertainty about its 
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true distance. We consider the possibility that the HI 
recession velocity may not give an accurate distance for 
this system, which could affect the derived absolute prop¬ 
erties of the objects in the HI1232+20 system. 

Large scale peculiar motions have been observed 
around the Virgo Cluster of galaxies, even beyond its 
virial radius, due to its si gnificant gravit ational influ¬ 
ences. Recently, Karachentsev et al. (201 4j) used a l arge 


sample of 1801 galaxies (mostly from Karachentsev & 
Nasonova (2010) and with some new observations) in 


the vicinity ot the Virgo Cluster which have indepen¬ 
dent distance measurements (e.g., Tully-Fisher, TRGB, 
Cepheid) as well as measured recession velocities. This 
kinematic sample of galaxies was used to map out the 
zero-velocity surface around the cluster, which encloses 
the region of s pace where gala xies are fa lling into the 
Virgo Cluster. Karachentsev et al. (20141 find that the 
zero-velocity surface radius is 7.2±0.7 Mpc, which corre¬ 
sponds to a projected radius of 25° ±2° at their assumed 
Virgo distance of 11=17.0 Mpc. The HI1232+20 system 
is at a projected distance of only ~8° from the center 
of the Virgo Cluster (NGC 4486), and may be partici¬ 
pating in t he inf all motion. Figure 1 in |Karachentsev fe| 
Nasonova (20101 shows a graphical representation of the 
difficulty in determining distances from recession veloc¬ 
ities in this region, and for this system’s measured re¬ 
cession velocity of ~1300 km/s, there are three possible 
distances. If HI1232+20 were infalling from the near side, 
located at the center of the cluster, or infalling from the 
far side, it could have distances of ~12, ~17, or ~25 Mpc, 
respectively. Further complicating the nearby velocity 
field is the Coma I cloud just north of the HI1232+20 
system. This complex of galaxies with peculiar velo cities 
is centered around (a, S) = (12.5h, -(-30°) (Karachentsev 


et al. 2011). Still, we can make a crude but reliable es¬ 


timate ot the lower limit on the distance based on the 
fact that we do not resolve any individual stars in the 
optical counterpart of AGC 229385. WIYN has been 
used to successfully resolve stellar populations in galax¬ 
ies under similar obser vin g con ditions out to 2—4 Mpc 
(e.g., Leo P, Rhode et al. 2013 M81 group, Rhode, pri¬ 
vate communication). Furthermore, WIY N observations 
of SHIELD galaxies (Cannon et al. 2011) resolve upper 
main sequence and supergiant stars at distance of 8 Mpc. 

Many of the most extreme properties of the objects in 
the HI1232+20 system are distance-independent quan¬ 
tities and would not be affected by a more nearby dis¬ 
tance (e.g., M hi /L b , surface brightness measurements 
and limits, average HI surface density). However, if they 
were only ~12 Mpc distant instead of 25 Mpc, the ab¬ 
solute quantities (e.g., HI mass, stellar mass, total lumi¬ 
nosity, physical area) would all scale down by a factor of 
four. For example, the HI mass of AGC 229385 would 
become 1.8 x 10 8 M e , its stellar mass would be reduced 
to 3.8 x 10 5 Mq, its absolute B magnitude would be re¬ 
duced to M b =— 11-4 mag, and its optical diameter would 
be ~2 kpc. 

Even with the difficulties of constraining the abso¬ 
lute distance to the HI1232+20 system, we are inter¬ 
ested in the large scale environment around it, and how 
isolated it has been on longer time scales. The back¬ 
ground spiral galaxy AGC 222741 (CGCG 129-006, la¬ 
beled in Figure!!]) appears on the sky between the three 
sources in HI1232+20, but has a recession velocity of 


i>i sr =1884 km/s, which is substantially higher than the 
velocities of the three sources in the system (1277, 1309, 
and 1348 km/s), so is very unlikely to be related. Fig¬ 
ure [6] shows the galaxies in the area around these HI 
sources. Galaxies shown on the plot come fro m the Up¬ 
dated Zwicky Catalog (UZC,|Falco et al.|1999|), from th e 
spectroscopic sample of SDSS DR9 (|AJin et al.||2012|), 
and from the Arecibo General Catalog (AGC). The UZC 
is an extragalactic redshift survey of ~20, 000 galaxies 
that is 96% complete to mz w <= 15.5 mag and the SDSS 
DR9 spectroscopic survey includes spectra of ~1.5 mil¬ 
lion galaxies, and is 95% complete down to r'= 17.7 mag. 
On Figure [6] a red “x” indicates the location of the back¬ 
ground spiral galaxy (AGC 222741), which appears near 
our HI sources on the sky. Also shown on Figure [6] is 
the projected virial radius of the northern sub clus ter of 
the Virgo Cluster of galaxies (Binggeli et ah 19851). The 


curves enclose a region centered at the position" of M87 
with a radius of 5.4° and centered on a recession veloc- 
ity of 1100 km/s and extending 800 km/s on either side 
( Ferrarese et al.|[2012 |. 

In order to demonstrate the complexity of the veloc¬ 
ity field around HI1232+20, we calculate 3-space separa¬ 
tions between this system and any nearby sources, using 
only their positions and observed velocities. The three 
sources within the HI1232+20 system are all within ~20' 
and ~70 km/s of each other. While the global veloc¬ 
ity field around these sources is compli cated due to t he 
influ ence of t he nearb y Virgo Clu ster ([Ka rachentsev fc 
Nasonova [2010 Karachentsev et al. 201 4[). we use this 
approach to crudely identify any possibfegalaxies which 
may be near enough to affect this system, and have sim¬ 
ilar positions and velocities. The object with the most 
similar velocity and position is NGC 4561 with a ve¬ 
locity of 1360 km/s and an angular separation of 1.4°, 
implying a physical separation of ~400 kpc assuming 
a simple Hubble flow. However, the Tully-Fisher dis - 
tance to NGC 4561 is 12.3 Mpc (Tully & Fisher 1987), 
so it is likely infalling to the Virgo Cluster from the near 
side. The object with the next-smallest 3-space distance 
is the starbursting galaxy IC 3605, which has a velocity 
of 1360 km/s and is located ~1.8° to the SE, with an im¬ 
plied physical separation (Hubble flow only) of ~500 kpc. 
No other distance measurements exist for IC 3605. The 
locations of NGC 4561 and IC 3605 are indicated on Fig¬ 
ure [6] with large dotted circles. The close proximity of 
this system to the Virgo Cluster means that its distance 
is uncertain and that the global velocity field is rich and 
complex. 

4. DISCUSSION 

The objects in the HI1232+20 system are not easily ex¬ 
plained, and some of their properties seem contradictory 
and puzzling. For example, it is difficult to understand 
the star formation history of AGC 229385, which appar¬ 
ently has only produced a tiny population of stars in an 
otherwise massive HI cloud. The HI mass of AGC 229385 
(log Mhi = 8.9) is larger than the nea rby star-f o rming 
Large Magellanic Cloud (log Mhi = 8.7,|Kim et al. 1998|) 
and just sm aller than the nearb y spiral galaxy 1V133 (log 
Mhi = 9.1, Gratier et al.|20 10), although its stellar pop¬ 
ulations and optical luminosity are vastly dissimilar. The 
HI kinematics of its HI cloud are also perplexing, as its 
rotation speed seems inadequate for its large mass and 
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Figure 7. Relationship between the HI mass-to-light ratio 
[Mhi/Lb, in solar units) and absolute B magnitude. The small 
dots and contours indicate Hi-detected galaxies from the AL¬ 
FALFA a.40 catalog matched with SDSS DR7 photometry. The 
dark circles and squares come from the new and archival obser¬ 
vations of|Warren et al.| (|2007|). The green points show the five 
(almost) dark galaxies from the pilot observations of |Cannon et al. | 
(2015). AGO 229385 is shown as a red star, and the upper lim- 
its tor AGC 229383 and AGC 229384 are indicated with diamond 
points and arrows. The error bars on AGC 229385 also include 
a distance u ncertainty of ±5 Mpc. The object from the |Ca.nnoli] 
|et al. ( 2015) sample with a larger lower limit on Mm/L b than 
the measured value for AGC 229385 is AGC 208602 and is likely a 
tidal feature and not an isolated galaxy. 


size. It is similarly difficult to explain why the other 
members of the HI1232+20 system have not formed any 
detectable stars, even with their substantial, although 
quite spread out, HI distributions. For comparison, both 
AGC 229383 and AGC 229384 have larger HI masses 
than the n earby dl galaxies IC 10 and NGC 6822 (lo g 


M hi = 8.0, Nidever et al.12013 


de Blok & Walter| 2000), 


but lack any optical counterparts m our observations. 

In order to put these objects in context with other 
galaxies, we consider their locations in typical galaxy 
scaling relations. In this exercise, we are treating these 
objects as independent galaxies and not simply gas 
clouds that have been stripped or tidally disturbed. It 
is possible that objects like this may be part of a large 
but mostly -unobserved cla ss of galaxies (e.g., the sunken 
galaxies of [Disney (1976|), but is more likely that the 
objects in this system are simply unique and uncommon 
galaxies. Since AGC 229385 has an optical counterpart 
we consider its location on optical and HI scaling rela¬ 
tions, while for the other sources our optical upper limits 
can still help constrain some of the same scaling relations. 
After discussing these scaling relations, we will will con¬ 
sider some possible formation scenarios to explain this 
unusual system of objects. 


4.1. Mhi/L relationship 

One of the most extreme properties of the objects 
in the HI1232+20 system is their exceptionally large 
HI mass-to-light ratio measurements (or lower limits). 
Galaxies are known to follow a typical relationship be¬ 
tween between this HI mass-to-light ratio and the over¬ 


all luminosity. The HI mass-to-light ratio is defined as 
Mhi/L , where L is the optical luminosity, often mea¬ 
sured in a blue filter. This relationship is especially dif¬ 
ficult to measure for faint low-mass galaxies, where a 
significant fraction of the optical luminosity may come 
from low surface brightness regions. Almost universally, 
whenever a galaxy with a reportedly large Mhi /Lb ratio 
is observed with deeper optical imag es, t he ratio returns 
to more typical values near unity. [Warren et al. ([20041 
used survey and catalog data to identify possible galax¬ 
ies with a large M H i/L b ratio {3<M H i / L B <27), but 
their sample of 9 large M B i/L B galaxies were almost 
all found to have less extreme r atios (Mht/L b < 5) after 
deeper observations. Similarly, (van Zee et al. 1997) used 
broadband optical imaging observations of six low surface 
brightness dwarf galaxies to show that their catalogued 
optical magnitudes had been severely underestimated by 
~1.5 mag, so their previously reported M B i/L b ratios 
became 4 times smaller and less extreme. Among dwarf 
galaxies, typical measurements of M B i/L b are between 


Sm/Im galaxies (Roberts & Haynes 

1994.; 

Stil & Israel 

2002) and held dls (Lee et al. 2003). 

The re 

iativeiy small 


est and highest ratio) highlights the importance of careful 
and accurate measurements of this ratio. 

Figure [7] shows the relationship between the HI mass- 
to-light ratio ( Mhi/L b ) and absolute B magnitude 
{M b ) for multiple samples of galaxies. The main dataset 
shown in Figure [7] comes from t he mat ched ALFALFA 
a.40 and SDSS catalogs of|Haynes et al.| (2011), which is 
the parent sample from which the ALFALFA (Almost) 
Dark Galaxies are drawn. We note that the shallow 
SDSS photometry may underestimate the luminosity of 
the faint sources in this sample. The average exposure 
time in SDSS is only ~1 minute, and will not detect 
low surface brightness emission from these galaxies. Ad¬ 
ditionally, optical fluxes may be erroneously estimated 
due to issues with SDSS background subtraction and ef¬ 
fects from nearb y bright stars. Also shown on Figure [7] 
are galaxies from |Warren et al. (2007), who observed 38 
galaxies at optical and radio wavelengths, and also com¬ 
piled an archival sample of previously observed galax¬ 
ies. Some of these sources are observed at low Galactic 
latitude, which results in additional uncertainty in the 
necessary extinction corrections. The five (al most) dar k 
galaxies f rom the pilot VLA observations of |Cannon et| 
al. (2015) are also shown on Figure [7] as green dots. The 


VLA observations allowed OCs tone identified for all 
but one of their sample, and the gas mass to light ra¬ 
tios for these sources are well-measured. We find that 
AGC 22 9385 has Mrrr/Ln /=45.8M^/L^, or, converted 
to B via|Jester et al.|(2005l, Mhi /Lb=38.2 Mq/Lq, and 
its position is indicated on Figure [7] As AGC 229383 and 
AGC 229384 are not detected in our deep optical images, 
we can only determine lower limits on Mhi / L g i , and find 
>31 and >57 Mq/Lq , respectively. These limits are also 
shown in Figure [71 where we have converted our upper 
limits in M g f to M B assuming the same g' — r' color as 
AGC 229385. 


Warren et al. (2007) suggest that there is an upper en¬ 
velope in Figure [7J which may represent the minimum 
amount of stars a galaxy will form, given a shallow po¬ 
tential well and an isolated environment. The sources 
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in the HI1232+20 system are in an extreme region of 
Figure [ 7 ] near this upper envelope. While there are 
other galaxies from ALFALFA with more extreme val¬ 
ues of Mhi/Lb shown on the plot, none have ratios 
that are as well-determined as the objects in this sys¬ 
tem. SDSS photometry for faint low surface brightness 
galaxies will likely underestimate their luminosity, which 
leads to an overestimate of the Mhi/Lb ratio. A recent 
study identified a low surface brightness galaxy near the 
Virgo Cluster with a very large HI mass-to- light ratio 


that they m easure as Mni/Ly > 20M 0 /L 0 (Bellazzini 
>jb AGC 229385 has the largest accurately 


et al. 2015 

measured HI mass-to-light ratio in the literature, but 
still appears to lie along a con tinuation o f the tren d seen 
in more luminous galaxies. Warren et al. (2007) use a 
similar sample of galaxies with Hi and optical data to 
fit the relationship between Mhi /Lb and Mb ■ They fit 
the upper envelope of the relationship with the following 
expression: 

log (M hi /L b ) max — 0 * 19(M b + 20.4). 

For AGC 229385 its Mb=— 12.72 would predict a maxi¬ 
mum Mhi/Lb= 29. While AGC 229385 does follow the 
general trend of low stellar mass objects having higher 
gas mass-to-light ratios, we measure M H i/L b = 38, w hich 
is even more extreme than the upper envelope of Warren] 
eFah]((2007). 


4.2. Galaxy Scaling Relations with Stellar Mass 

Studies of large samples of galaxies have found that 
stellar mass seems to be an important parameter that re¬ 
lates to star formation in the possible evolution of galax¬ 


ies from the blue cloud to the red sequence (Brin chmann 
et al.|2004 Salim et al. 12007 Huang et al.|2012b). Using 
a sample of 9,417 ALFALFA-selected galaxies with coun- 


sample of 9,417 ALFALFA-selected galaxies ■ _ 

ter parts in archival GALEX and SDSS images, |Huang et| 
al. ( 2012 b) studied the scaling relations as a function of 
stellar mass and optical color (see their Figure 8 ). Fig¬ 
ures and [ 8 j; show the HI mass and fni—M hi/ 'M* as 
a function ofthe stellar mass, which is determined from 
SED fits. A clear relationship with Mhi is found from 
M*=3.2 x10 7 Mq to 3.2 x10 11 Mq, with a change in slope 
at M*~ 10 9 M q . Analogously, f hi follows the same gen¬ 
eral trend with a break at ~1O 9 M 0 . On both panels the 
location of AGC 229385 is indicated with a large star, 
and the upper limit measurements of AGC 229383 and 
AGC 229384 are indicated with arrows. These sources 
are deviant from the expected scaling relations at the low 
stellar mass end, and have too much HI for their stellar 
mass (detected or not). 

Figures [ 8 ]:) and[ 8 ]i show Mhi and fni as a function of 
NUV — r color, which acts as an indicator of the amount 
of recent star formation (UV) compared to the amount 
of past star formation (r). Here only AGC 229385 can be 
plotted, since we have no optical or UV detections of the 
other two sources. AGC 229385 is at an extreme location 
in both of these parameter spaces, and lies significantly 
above an extrapolation of the low mass trend in the rela¬ 
tionship between fni and NUV—r. The color cannot be 
much bluer than it already is, since after 5 Myr, a sim¬ 
ple stellar population_of half solar metallicity will have 
NUV— r=—0.2 (Bressan et al. 2012). AGC 229385 seems 
to be deviant in the sense thafTt has too much HI for 
its NUV — r color, consistent wit h the prev iously dis¬ 
cussed panels. Finally, we note that Toribio et al. (2011) 


a.40-SDSS-GALEX 



Figure 8. The scaling relations between HI and optical obser¬ 
vations, using the O.40-SDS S-G ALEX sample and analysis from 
Figure 8 of (Huang et al.| (|2012b|. In all panels, number densities 
of galaxies are calculated within grid cells set by the intervals of 
the minor ticks on the axes. Contours are drawn above densities 
of 20 galaxies per grid cell, and selected contours are labeled with 
their number densities. Blue diamonds and lines show the average 
y-values in each x-bin. Panel (a) shows the relationship between 
HI mass and stellar mass, with stellar masses from mass-to-light 
ratios from fits to SDSS photometry. Panel (b) shows the HI mass 
as a function of NUV-r color. Panel (c) shows the ratio of HI mass 
to stellar mass Panel (d) shows the relation¬ 

ship between fui to NUV-r color. Typical error bars of individual 
galaxies are shown in the corner of panels (a) and (d). The red star 
shows the location of AGC 229385 in each panel, while the red ar¬ 
rows indicate the upper limits of AGC 229384 and AGC 229383 
where measured. 

found a strong relationship between Mhi and the opti¬ 
cal diameter measured at p = 25mag/arcsec 2 , such that 
galaxies with smaller values £>25 typically had smaller HI 
masses as well. However, we cannot compare the prop¬ 
erties of the HI1232+20 system with this relationship 
because the optical counterpart to AGC 229385 never 
reaches /i = 25mag/arcsec 2 , and the other two sources 
have even fainter limits on their optical non-detections. 


4.3. Galaxy Scaling Relations with HI mass 

We now discuss the scaling relations which depend on 
the total HI mass of a galaxy. First, we consider the rela¬ 
tionship between HI diameter (measured at a parti cular 
t he co lumn density) and total HI mass. Broeils & Rhee 
( 1997| hereafter BR97) observed this relationship in a 
sample of 108 spiral and irregular galaxies and found a 
strong correlation between logMni and logDni (the HI 
diameter at lM 0 /pc 2 ), with a d ispersion of o nly 0.13 
dex, and a slope of 1.96 ± 0.04. |Swaters et al.|(2002| 
found the same correlation and slope when they used a 
sample of lower mass irregular galaxies. This relation¬ 
ship has been measured between HI masses of 6 x 10 7 
and 3 x 1O 1 O M 0 and between HI diameters of 0.8 and 
160 kpc. A consistent relationship implies that there is 
a constant average HI density in all gas-rich galaxies, 
which BR97 estimates as 3.8 ± l.lM 0 pc -2 . 

While explanations for the underlying mechanisms for 
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this correlation are complex (e.g., feedback, turbulence, 
etc.), it is simple to apply the relationship from BR97 
to our sources. In the case of AGC 229385, our HI 
mass predicts an HI diameter of 16 kpc (measured at 
lAf 0 /pc 2 ). Our HI observations of AGC 229385 show 
an elongated source with dimensions 24 x 10 kpc, when 
measured at this HI column density. This is ~50% larger 
than predicted, and while projection effects could re¬ 
duce this discrepancy somewhat, the HI kinematics do 
not suggest we are viewing a disk edge-on so any cor¬ 
rection for projection effects will be small. AGC 229384 
has an HI mass of 2.00 x 1O 8 M 0 which predicts an HI 
diameter of 8.1 kpc. Our observations show an elon¬ 
gated HI source with dimensions 10 x 8 kpc, which is 
slightly larger than expected (although this source has 
two separate HI peaks at that column density). The 
HI distribution of AGC 229383 never reaches the column 
density that corresponds to lM 0 /pc 2 , and we are unable 
to apply this relationship to it. Both AGC 229385 and 
AGC 229384 have larger HI diameters than predicted by 
the HI diameter-HI mass scaling relation. This also im¬ 
plies that their average HI surface density is significantly 
less than the constant value found by BR97. Assuming 
a circular HI disk for AGC 229385 and AGC 229384, we 
find that their average HI surface densities are 1.6 and 
2.4M 0 pc -2 , respectively, which ar e at the lo wer limit of 
the di stribution found by BR97. Rosenberg fe Schneider| 
(2003) also compared Mjji and Djji for a different sam¬ 
ple” but their relationship gives similar results to BR97. 
These unusually low HI surface densities may be related 
to the lack of significant star formation in the objects in 
the HI1232+20 system. 

4.4. Galaxy Scaling Relations with HI kinematics 

We consider the HI kinematics of the sources in the 
HI1232+20 system, especially with regard to their ap¬ 
parently slow rotation. However, since the HI rotation 
curves of these sources are difficult to fit or interpret (see 
Figures 3^, Hi, andj5^), we instead use the integrated 
width o ft lie 21 cm line itself to measure their rotation. 
The HI velocity widths (W50, measured at 50% of the 
peak value) of the three sources in the HI1232+20 system 
are unusually small for their HI masses. We model the 
HI velocity width distribution for all ALFALFA galax¬ 
ies as a function of their HI mass. Integrating over 
the model distribution at the HI mass of AGC 229385 
(Mhi= 6.7 x 1O 8 M 0 ) we find that only 2% of ALFALFA 
galaxies have velocity widths smaller than its measured 
value of W 5 o =34 km/s. Similarly for AGC 229384, we 
find only 3% of ALFALFA galaxies at its HI mass have 
similarly small values of W 50 . We note that the much 
wider velocity width of AGC 229383 falls near the 50th 
percentile of objects of its HI mass, but that this may be 
due to the presence of multiple objects within the AL¬ 
FALFA beam. Full details of the velocity width model 
will be published in a forthcoming paper (Jones et al. in 
prep.). 

We can next use the Baryonic Tully Fisher rela¬ 
tion (McGaugh 2012) to compare the total baryonic 
masses with the rotation velocities for the members 
of the HI1232+20 system. Since these o bjects are 
gas-d ominated, we use the relationship from |McGaugh] 
( 2012 ) to calculate rotation velocities as 
where W 20 is the HI velocity width at 20% of the maxi- 



Figure 9. Bary onic Tullv- Fishcr relation. Black do ts sho w galax¬ 
ies measured by |McGaugh| (|2012|l and |McGaugh| (|2005f. Open 
polygons show the sources in the HT12.12+211 system from this 
work. AGC 229385 (shown as an open pentagon) has the largest 
HI mass of this sample and is offset above the BTF relation. The 
next most massive, AGC 229384 (shown as an open square), is also 
offset above the relationship, while the least massive, AGC 229383 
(shown as an open triangle), is consistent with the BTF relation. 
The two HI peaks of AGC 229383 are measured separately (shown 
as small triangles), and both lie along the BTF relation as well. 
The “x” shows the position of Leo P, a metal-poor gas-rich dwarf 
galaxy discove red by ALFALFA in the local volume just outside the 
Local Group llGiov anelli et al. 12013 Rhode et al.|2013| |Bernstein-| 
|Cooper et a l. |zU14[ l. ' 

mum. While the HI rotation of these objects is difficult 
to measure accurately, their placement on the BTF re¬ 
lation is intended as a suggestive exercise to shed light 
on some of their unusual properties. If these objects are 
indeed galaxies, then their rotation seems too slow for 
their measured mass. 

Figure [9] shows the locations of the objects in the 
HI1 232+20 system compared to a la rge sample of galax¬ 
ies ( McGaugh|2012 McGaugh 2005). The baryonic mass 
is the sum of the stellar and total gas mass, and in gas- 
dominated galaxies without substantial stellar popula¬ 
tions, is determined as Mb ary =1.33 x Mri, where the ex¬ 
tra factor of 1.33 accounts for helium. The stellar mass of 
the optical counterpart of AGC 229385 contributes only 
a negligible ~0.2% compared to its HI mass. The upper 
limits on the non-detections of stellar populations of the 
other two objects indicate that these would contribute 
similarly negligible amounts of baryonic mass. This es¬ 
timate of baryonic mass does not include contributions 
from molecular or ionized gas. It is conceivable that there 
may be an envelope of low density ionized gas surround¬ 
ing these galaxies, contributing more mass than is in¬ 
cluded in our determination of Mb ary . 

Figure [9] shows that the two more massive objects 
in the HI1232+20 system both fall significantly above 
the standard BTF relationship, while the least mas¬ 
sive (AGC 229383) is in agreement. However, since 
AGC 229383 has two strong HI peaks which are sep¬ 
arated by ~ 40 kpc, we also consider the kinematics of 
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each clump separately, and plot them on Figure[9]as well. 
AGC 229383 is the lowest column density object of this 
system, and is, at best, difficult to interpret. Naive place¬ 
ment of the combined object or its peaks seems to agree 
reasonably well with the BTF, but the low signal-to- 
noise nature of the HI observations makes further anal¬ 
ysis or interpretation difficult. We also caution that the 
WSRT observations do not recover the total HI flux of 


ALFALFA (see Section 3.3), so the HI masses of the two 
peaks do not sum to the total observed mass. 

We note that inclination effects could slightly mod¬ 
ify our determination of rotation velocities, but that the 
discrepancies from the BTF are larger than can be ac¬ 
counted for by changing the inclination angle of the ob¬ 
jects. If the distance to the HI1232+20 system was signif¬ 
icantly smaller, then they would be in better agreement 
with the BTF. In order for AGC 229385 to agree with the 
BTF it would need to be at a distance of ^4.4 Mpc, and 
AGC 229384 would need to be at a distance of ~5.4 Mpc. 
This smaller distance is unlikely given other constraints 
(e.g., the lack of resolv ed stars in our optical observa¬ 


tions, see Section 3.4.2). 

We also consider the ratio of dynamical mass to 
HI mass (M dyn /M m ) as another independent con¬ 
straint on the distance by assuming a typical value 
of Md yn /AlHi=10. We measure an effective velocity 
(VeS~ v iot + 3cr 2 ) for AGC 229385 of v e s=2\.2 km/s, so 
this relationship implies a distance of 4.3 Mpc. If we 
assume no dark matter (e.g., M dyn /MHi=1.3), we find 
a distance of 32.8 Mpc instead. The HI kinematics of 
AGC 229384 and AGC 229383 make it difficult to mea¬ 
sure v ro t, which is required to determine v e s. 

4.5. Formation scenarios 

Given the variety of observational constraints we have 
for the sources in the HI1232+20 system, and the context 
from existing galaxy scaling relations, we now consider 
some of the possible evolutionary scenarios which could 
account for a system like this. 

The blue color of the optical counterpart of 
AGC 229385 suggests that it may consist of a mostly 
young stellar population. If this object has only just 
begun forming stars, it may not have had enough time 
yet to convert a significant amount of its gas into stars, 
as reflected by the high gas mass-to-light ratio. How¬ 
ever, the lack of Ha detection and the weak UV SFR 
(~O.OO4M 0 /yr) seem to indicate that it is only slowly 
forming stars. At this rate it would take more than a 
Hubble time to generate enough mass in stars to return 
it to the normal relationship between HI and stellar mass 
(e.g.. Figure^). 

It is difficult to find a single convincing explanation for 
why one HI cloud has an observable stellar population 
while its two nearby neighbors do not. Obviously, this 
creates a “fine-tuning” problem of sorts. Interactions 
between the members of this system may have triggered 
star formation in the most massive object. It is also 
possible that the system may have been pertur bed by 
an external object. Our analysis in Section [3.4. 1| showed 
that there was only one plausible perturber nearby, even 
given our generous assumptions about relative velocities 
and timescales for interactions. 


Alternatively, Verde et al. (2002) has suggested that 
low-mass dark matter halos can contain neutral gas with¬ 


out ever forming stars, under certain conditions. How¬ 
ever, the members of this system have substantial HI 
masses and are likel y differ ent from the low-mass dark 
matter halos of Verde et al. (2002). It may be that 
star formation has only occurred in the HI cloud which 
is dense enough at its cen ter to exceed a gas density 
threshold (e.g., Kennicutt (1998)). Indeed, in studies 
of the outer disks of larger galaxies, star formation ap¬ 
pears to cease below an HI surface density of ~ 1M 0 /pc 2 
(Radburn-Smith et al.|2012 Hunter et ai.|2013 1. The HI 
distributions of the sources m the H1I232+20 system are 
mostly below this surface density threshold, and only the 
peak of AGC 229385 substantially exceeds it. 

The reasons behind the low gas surface density and 
over-extended HI distributions of these objects are not 
clear. In the case of AGC 229385 it is possible that feed¬ 
back from the star formation (perhaps stronger in the 
past) has injected kinetic energy into the HI, and tem¬ 
porarily expanded it. Or, we may be seeing the results 
of a recent infall of cold gas into this system. However, 
these objects have significant amounts of HI, and it is 
difficult to justify an inflow scenario which could pro¬ 
vide enough gas while maintaining a low gas density and 
inhibit star formation. 

It is also possible that these objects may be the most 
massive knots in a system of smaller HI clouds which 
are in the midst of a tidal interaction or merger event. 
AGC 229385 and AGC 229384 may be the brightest HI 
peaks, and could be the products of recent mergers of 
smaller HI clumps. The extended nature of AGC 229383 
and its two HI peaks are suggestive of an extended HI dis¬ 
tribution with local peaks and clumps. However, the AL¬ 
FALFA observations are very sensitive to faint HI emis¬ 
sion, and given the good agreement between WSRT and 
ALFALFA HI fluxes for AGC 229385 and AGC 229384, 
it seems unlikely that there is a substantial amount of 
unseen extended gas in most of the system. 

The HI1232+20 system does not seem to be a tidal 
feature or remnant of a larger object, based on the lack 
of a clear connection to any possible external perturb¬ 
ing object. However, tidal interactions between mem¬ 
bers of the system may be important to their individual 
star formation histories. The HI clouds in this system 
have extremely low star formation efficiencies and low 
gas densities. These low gas densities are especially dif¬ 
ficult to reconcile with the significant HI masses of these 
sources. The HI mass of AGC 229385 is greater than 
that of the LMC, and AGC 229383 has more HI than 
some nearby star-forming dwarf irregular galaxies. Even 
if recent tidal interactions between the sources may have 
triggered a burst of star formation in AGC 229385, the 
overall properties of the HI1232+20 system are difficult 
to convincingly interpret and explain. 


5. SUMMARY 

In this work we present the discovery of the HI1232+20 
system of HI sources, drawn from the sample of (almost) 
dark extragalactic sources in ALFALFA. This system 
defies conventional explanations and our HI synthesis 
imaging and deep optical observations have revealed a 
set of objects with properties that are difficult to rec¬ 
oncile with typical scaling relations. The most mas¬ 
sive of its members (AGC 229385, Mhi= 7.2x 1O 8 M 0 ) 
has a weak stellar counterpart, detected in UV and 
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ultra-deep optical imaging, with a peak surface bright¬ 
ness of /r s /=26.4 mag/arcsec 2 . It has the most extreme 
well-measured gas mass-to-light ratio in the literature 
(. Mhi/Lb = 38), and its absolute magnitude is only 
Mg* = —12.9 mag (M* = 1.5x10 6 Mq), assuming a 
distance of 25 Mpc. We do not detect optical counter¬ 
parts for the other two members, but place upper lim¬ 
its on their absolute magnitudes of M g > > —11.3 mag 
( M+ < 3x 10 5 Mq). The HI kinematics of the three ob¬ 
jects in this system are inconsistent with typical galaxy 
scaling relations, with HI distributions that are too ex¬ 
tended and too slowly rotating for their HI mass. This 
group appears on the sky just outside of the projected 
virial radius of the Virgo Cluster, but is otherwise iso¬ 
lated from any nearby galaxies. 

The HI1232+20 system is difficult to explain com¬ 
pletely, but may be an example of objects just above and 
just below a threshold for star formation. The most mas¬ 
sive of the three sources is forming stars, but may have 
only recently started to do so. The other two sources 
have no observational signatures of star formation, so 
there may be some mechanism inhibiting this process. 
Sources like these are very rare in the ALFALFA survey, 
especially at such large HI masses. As observations of 
the HI1232+20 system continue we hope to learn more 
about its history and role in galaxy formation and evo¬ 
lution. 
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